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We propose a method to achieve coherent couphng between Nitrogen-vacancy (NV) centers in 
diamond and superconducting (SC) flux qubits. The resulting coupling can be used to create a 
coherent interaction between the spin states of distant NV centers mediated by the flux qubit. Fur- 
thermore, the magnetic coupling can be used to achieve a coherent transfer of quantum information 
between the flux qubit and an ensemble of NV centers. This enables a long-term memory for a SC 
quantum processor and possibly an interface between SC qubits and light. 
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Among the many different approaches to quantum in- 
formation processing, each has its own distinct advan- 
tages. For instance, atomic systems [l[ present excellent 
isolation from their environment, and can be interfaced 
with optical photons for quantum communication. In 
contrast, condensed matter systems [1] offer strong in- 
teractions, and may benefit from the stability, robust- 
ness, and scalability associated with modern solid-state 
engineering. In the last years, much effort is being de- 
voted to coupling atomic and solid-state qubits to form 
hybrid systems, combining "the best of two worlds". One 
attractive approach to hybrid systems involves the inte- 
gration of atomic ensembles with superconducting (SC) 
stripline resonators. Strong coupling between SC qubits 
and such resonators has already been achieved and 
approaches to extend the coupling to atomic systems 
have been proposed [1,01 • To achieve an appreciable cou- 
pling, these proposals often use an electric interaction, 
but magnetic interactions are more desirable, since long 
coherence times are mainly achieved in systems where 
spin states are used to store the information. Magnetic 
interactions are, however, inherently weaker but recently 
it has been proposed theoretically 0] and shown ex- 
perimentally 0] that strong coupling to ensembles of spin 
systems can be achieved. 

In this Letter, we propose a novel hybrid system, 
consisting of a SC flux qubit magnetically coupled to 
Nitrogen-vacancy centers (NVs) in diamond. The latter 
system shares many of the desirable properties of atoms, 
such as extremely long coherence times and narrow-band 
optical transitions 0], but at the same time, the integra- 
tion with solid state systems can be relatively easy, as it 
eliminates the need for complicated trapping procedures. 
Additionally, much higher densities can be achieved with 
very limited decoherence @, Q . As we show below, the 
magnetic coupling between a SC flux qubit and a sin- 
gle NV center can be about three orders of magnitude 
stronger than that associated with stripline resonators. 
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FIG. 1: (a) Schematic setup. An isotopically pure ^^C dia- 
mond crystal, doped with Nitrogen-vacancy color centers, is 
located in the proximity of a flux qubit of size Lx L and cross 
section h x h. The two persistent-current quantum states of 
the flux qubit have different associated magnetic fields, which 
give rise to a state-dependent interaction with the electron 
spin in the NV center(s). (b) The combined system of a flux 
qubit and a NV center. The eigenstates of the flux qubit are 
superpositions of left- and right-circulating currents. An ex- 
ternal magnetic field splits the ms = ±1 states of the NV, 
resulting in a two-level system with the ms = ms = 1 
transition close to resonance with the flux qubit. 



thereby making the system an attractive building block 
for quantum information processing. 

Flux qubits (FQs) form superpositions of persistent 
currents of hundreds of nano- Amperes, flowing clock- 
wise and anti-clockwisc through micromcter-sizcd SC 
loops [l^. The magnetic field associated with this cur- 
rent, of the order of a ^T, enables a magnetic dipole 
coupling to the electron spins associated with crystalline 
impurities such as the NV center in diamond. Of particu- 
lar interest is the coincidence of energy splittings: the two 
states of the FQ are typically separated by a few GHz, 
while NV centers have a S" = 1 ground state, with zero- 
field splitting A = 27r X 2.87 GHz between the ms = 
and ms = ±1 states. By the application of a mT mag- 
netic field, one of the spin transitions of the NV center 
can be tuned into resonance with the FQ (see Fig. [T](b)). 
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This, together with the large magnetic moment of the 
FQ and the relatively long coherence times of both sys- 
tems, opens the possibility of achieving coherent transfer 
between them. 

Let us consider a single NV center in a diamond crys- 
tal, located near a square FQ of size L and thickness h 
(see Fig. [Ija)). A static external field B'^^'^ is applied, 
whose component perpendicular to the FQ provides half 
a flux quantum, and brings the qubit near the degener- 
acy of the clockwise and counter-clockwise current states. 
The zero-field spin splitting of the NV center A sets a pre- 
ferred axis of quantization to be along the axis between 
the Nitrogen and the vacancy; thus, the field parallel to 
this axis sets the small additional Zeeman splitting be- 
tween ms = ±1 states, and allows to isolate a two-level 
subsystem comprised by ms = 0, 1. 

For convenience, we denote as the z axis the crystalline 
axis of the NV center. The Pauli operators for the FQ 
system, not tied to a particular spatial axis, will be de- 
noted by Ti, f2, fa, with fa describing the population dif- 
ference between the two persistent-current states. The 
interaction of the total magnetic field B (external and 
from the FQ) with the NV center can be written as S-W ^ 
where W = geHsB, Qe is the electron g- factor and /xs is 
the Bohr magneton; the two persistent-current quantum 
states of the FQ give rise to different anti-aligned mag- 
netic fields: T^W^'~^. The Hamiltonian for the system is 
then 



H = ef3/2 + Afi + A^f + + hW^^ ■ S. (1) 

Here the magnitude of W^'^ corresponds to a Larmor fre- 
quency shift due to the circulating or counter-circulating 
currents in the FQ, A is the coupling between these two 
current states, and e is the bias in the two- well limit of the 
FQ, which depends on the external field perpendicular to 
the loop. If the flux qubit's plane is not perpendicular to 
the z axis of the NV center, both systems can be tuned 
on resonance by changing independently the z and, e.g., 
X components of i?"'^' . 

Due to the coupling A in equation ([T]), the eigenstates 
of the FQ Hamiltonian arc not left- and right- circulating 
current states. This means that there is a magnetic tran- 
sition between the dressed states of the FQ which couples 
to the electronic spin of the NV. To describe this, we ro- 
tate the FQ via a unitary transformation by an angle 
cos Q = ejlu}^ giving two FQ dressed states with a transi- 
tion frequency a; = yJe'^/A + A^. When A -I- W°^^ — lo^S 
is small, we can transform to a rotating frame and make 
the rotating-wave approximation (RWA) to describe the 
near- resonance interaction between the NV and the FQ. 
Neglecting the state ms = — 1, due to the external field 
moving it far out of resonance {\6\ ^ the cffec- 
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FIG. 2: (Color online). Coupling g between a FQ operated at 
the degeneracy point (e = 0) and a single NV center located at 
the position {x, y, z) in a reference frame with axes as in Fig. [1] 
but centered in tiie middle of the square FQ. The NV crystal 
axis is assumed to be parallel to one of the wires forming the 
flux qubit. The FQ has size L = l/im, thickness h — 60 nm, 
and a critical current of 0.5 /xA. 



five Hamiltonian describing the dynamics is 



RWA 



V2 



H.c. 



(2) 

where & are Pauli operators describing the NV m^ = 0,1 
electron-spin states, and t is in a rotated basis so that, 
e.g., fa describes the difference in the populations of the 
FQ dressed states. 

The coupling constant g = sin ewl'^/V2 depends on 
the field perpendicular to the NV axis, and it is maximal 
at the degeneracy point e = 0, where 9 = tt/2. Further- 
more, this point has a "sweet-spot" property: the ener- 
getics of the system is insensitive to small fluctuations 
of e, as the eigenvalues have zero derivative with respect 
to e at this point. This reduces the dephasing of the FQ 
due to stray magnetic fields and, e.g., paramagnetic spins 
in the diamond crystal [llj, and ensures that there will 
be no differential shifts of the resonance frequency of the 
NVs due to an inhomogeneous static field from the FQ. 
For the remainder of this Letter we will assume that we 
are working at the degeneracy point. Figure [2] shows the 
coupling g, along four different vertical lines. We eval- 
uate this using the magnetic field generated by a finite 
width square loop as given by the Biot-Savart law. For 
a FQ of size L = 1 /im and critical current 0.5 /lA, the 
coupling reaches g = 2tt x 12 kHz for a single NV center 
located at the center of the loop, which is about a factor 
of 1000 larger than the coupling achieved with stripline 
resonators |6|. This coupling g is however too small to 
achieve coherent transfer, since current T2 times in FQs 
are at best a few microseconds 
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We now show that the FQ can be used as a virtual 
intermediary, allowing to couple coherently two or more 
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NV centers with the same orientation and detuned 5 from 
the FQ. For large enough detunings, 5 » l/Tj^'^,^, we 
can adiabatically ehminate the excited state of the FQ, 
and the two coupled NVs have the states |11) and (|10) + 
\01))/^/2 shifted by an energy 2g^/6. In contrast, the 
states (|10)-|01))/y2and |00) are not shifted by the FQ. 
As a result, for a fixed time tx = 7r5/(4g^) an operation 
resembling ySWAP - an entangling operation - between 
two NV centers can be implemented. 

To analyze the gate fidelity, we notice that the coher- 
ence of FQs and NVs are subject to low-frequency noise. 
Larger coherence times arc typically obtained by remov- 
ing this noise during single-qubit evolution, using spin- 
echo sequences. Such sequences can be incorporated into 
the gate operation by following the prescription for com- 
posite pulses of the CORPSE family In particular, 
turning on and off the interaction with the FQ, e.g., by 
changing from a 7775 ~ +1,0 to ms — —1,0 superpo- 
sition, local TT pulses (described by (Tx) allow to realize 
a sequence Ex^/4_0/2<5'xO'x'^Ex0WO'^(T^Ex^/4_0/2: where 
(f> = 2 sin~^(l/-\/8) ~ 41.4°, Exg is an exchange-type in- 
teraction for a rotation 6, with 6* = tt a "SWAP" , and w 
is a wait for a time tw = ^x(l/2 — 20/7r), which makes 
the time between the tt pulses equal to twice the time on 
either end of the sequence. This approach is only sensi- 
tive to detuning errors in fourth order for both collective 
and individual noise on the two spins, effectively inte- 
grating a Carr-Purcell-type spin echo with the vSWAP 
operation. Furthermore, the total time tx spent inter- 
acting with the FQ during the sequence is equal to the 
prior case, thus inducing no additional overhead in the 
virtual coupling through the FQ. 

With a spin echo the coherence is often limited by 
energy relaxation, resulting in an exponential decay 



12| . whereas the NVs decay 



exp(-t/T2^^) for the FC 

as exp{-{t/T^^f) 0, [lj|. During the operation, the 
finite chance of exciting the detuned FQ leads to an in- 
duced decoherence rate 7 = 2(jr^/(5^Tj'^), resulting in a 
gate error ~ {{tx + 2tw)/T^^ f + jtx- Minimizing this 
expression we find an optimal detuning, which gives the 
optimized error probability - 2.2/{g^T^T^"^)^/^ [HI. 
For isotopically purified ^^C diamond [l5| with an opti- 
mistic T^^ « 20 ms, g « 27r X 12 kHz, and 5 ^s, 
the maximum achievable fidelity of the ~ ySWAP oper- 
ation is ^ 0.98 with an operating time tx + 2tw ~ 3.3 
ms. It is thus possible to achieve a high-fidelity coher- 
ent operation between two NVs separated by micrometer 
distances. This coupling may even be extended to NVs 
separated by large distances: If two FQs are strongly 
coupled (directly or through resonators) with a coupling 
exceeding the detuning S, two NVs residing in different 
FQs may be coupled through the dressed states of the two 
FQs. This results in a long-distance coupling between the 
NVs of roughly the same magnitude as derived above. 

Even though the FQ coherence times are much shorter 



than the coupling to a single NV center, it is possible 
to coherently transfer the quantum-state from the FQ 
to an ensemble of many {N) NVs by benefiting from a 
^/N enhancement in the coupling constant [l^ • Consider 
the diamond crystal depicted in Fig. [IJa), with density 
n of NV centers, each of them with a fixed quantiza- 
tion axis pointing along one of four possible crystallo- 
graphic directions. If the orientations are equally dis- 
tributed among the four possibilities, and the external 
field is homogeneous, a quarter of the centers can be 
made resonant with the FQ. The total coupling will then 
be given by the FQ - single NV interaction summed over 
the resonant subensemble, whose state is conveniently 
expressed in terms of the collective angular momentum 
operator J+ = (1/G) J2j 9j^+\ where the sum runs over 
the resonant NVs, and G = (^^ ISjP)^^^ is the collec- 
tive coupling constant. At the operating temperature of 
the FQ (tens of mK) the NVs are near full polarization, 
dz ~ (ctI"''') — —1, and the collective spin operators fulfill 
harmonic-oscillator commutation relations, thus having 
bosonic excitations in this limit. Tuning the system into 
resonance, the interaction between FQ and ensemble of 
NV centers is 



H.c. 



(3) 



whose dynamics can be complicated in general jlTl llSj , 
but close to full polarization takes place between the 
states |1)fq|0^)nv and |0)FQi+|0^)Nv, where |0^)nv 
corresponds to all NVs being in the ground state jlOj . 
One can thus reversibly transfer the quantum state be- 
tween the FQ and the collective excitations of the NVs. 

The collective coupling G as a function of the density 
n of NV centers, is shown in Fig. [3] for three different 
sized FQs. The coupling is obtained by summing the 
inhomogeneous coupling over NVs distributed in a crystal 
of size {2L)^. Taking the FQ coherence time to be ~ 5 /zs, 
we find that for a FQ with L = 5 /im, coherent transfer 
becomes possible at densities n > lO^^cm"'^, which have 
already been achieved in recent experiments [l^j. We 
note that it might be possible to increase the interaction 
strength with the FQ by designing circuits with higher 
critical currents. Furthermore, an important feature of 
the present approach is that strong coupling to the FQ 
is achieved with ensembles containing a relatively small 
number of spins in a few micrometers, which makes it 
easier to achieve fast and identical manipulation of all 
spins, e.g., in spin-echo approaches. 

The performance of the transfer of information be- 
tween FQ and collective ensemble will be limited by 
paramagnetic impurities present in the diamond crys- 
tal, which will interact via a dipolar coupling with the 
NVs encoding the collective quantum state. In the Sup- 
plementary Information we give a detailed discussion of 
decoherence induced by these impurities, which, due to 
the typically low Nitrogen to NV conversion efficiency, is 
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FIG. 3: (Color online). Coupling between a FQ and an 
ensemble of NV centers, as a function of the NV density n. 
Results are shown for three sizes L of the square flux qubit, 
and are obtained by summing the inhomogeneous coupling 
constant, calculated as in Fig. [2] over an ensemble of NVs 
located in a cubic diamond crystal of size Lc = 2L placed at 
a height h/2 below the FQ. In all cases, the NV crystal axis 
is assumed to be parallel to one of the wires forming the FQ, 
the width of the FQ is /i = 60 nm, and its critical current is 
0.5 



the nuclear spin state (m/) through a sequence of radio 
and microwave frequency pulses performing the evolution 
a|00) + /3|10) a\00) + /3|11) a\00) + ;3|01), where 
the states label the magnetic quantum numbers \msmi). 
This allows for a long-term memory in the system while 
other operations are performed, e.g., while the NVs in- 
teract with light for quantum communication. 

We have shown how to magnetically couple a supercon- 
ducting FQ to NV centers in diamond. This may be used 
to achieve strong coupling of distant centers or to trans- 
fer the state of a FQ to an ensemble of NVs. The latter 
could be used for long-term nuclear storage, and, using 
the strong optical transitions of the NVs, may enable an 
interface between superconducting qubits and light. 

We thank P. Forn-Dfaz, K. M0lmer, R. Schoelkopf, 
A. Imamoglu, K. Scmba and J. E. Mooij for discus- 
sions. This work was supported by NSF, CUA, DARPA 
and Packard Foundation, and grants FPU AP2005- 
0720, FTP 274-07-0080, FIS2009-08744 and CCG08- 
CSIC/MAT- 3775. At the time of submission a related 
preprint appeared [1^, which describes an NV coupled 
to a persistent-current loop integrated into coplanar res- 
onators. 



dominated by unpaired Nitrogen electrons in the sample 
lll |. For a FQ of size L = 5 /.jm and a density of n ~ 10^^ 



cm~'^ we estimate a coherence time of T2 =1.8 fis. This 
is sufficiently long for the infidelity induced by the para- 
magnetic impurities to be on the percent level, both for 
the transfer from the FQ to the spin wave and the nuclear 
storage discussed below. Furthermore the dccohercncc of 
the FQ induced by the impurities is negligible if we work 
close to the degeneracy point of the FQ, cos9 w 0. 

So far we have ignored the influence of the nuclear spin. 
The strong hyperfine interaction with lattice ^'^C will be 
detrimental to the presented schemes, but this can be 
overcome using isotopically purified ^^C diamond [isj . 
For the Nitrogen atoms forming the centers there are no 
stable isotopes without nuclear spin. This spin can, how- 
ever, be polarized by transferring the nuclear state to the 
electron spin, using a combination of radio-frequency and 
microwave pulses, followed by polarization of the elec- 
tron spin [2l| , or directly through optical pumping using 
excited-state couplings [22|. For ^'^N, the large quadrap- 
olar field (~ 5 MHz) from the NV center leads to quan- 
tization of the spin-1 nucleus along the NV axis. This 
suppresses spin-flip terms such that the hyperfine inter- 
action just acts as an additional parallel component to 
the external field. Due to the long (^ 1 s) relaxation 
time of the ^*N nuclear spin, this can be accounted for 
by an appropriate detuning of the magnetic field so that 
the nuclear spin docs not lead to dccohercncc. 

The nuclear spin can actually be turned into a valuable 
resource for long-term storage of quantum information. 
The electronic spin state (ms) can be transferred into 



[1] 
[2] 

[3] 
[4] 

[5] 
[6] 



[7] 



[9] 
[10] 

[11] 
[12] 



[13] 
[14] 
[15] 
[16] 

[17] 
[18] 
[19] 
[20] 



J. I. Cirac and P. ZoUer, Phys. Rev. Lett. 74, 4091 (1995). 

B. E. Kane, Nature 393, 133 (1998); D. Loss and D. 
DiVincenzo, Phys. Rev. A 57, 120 (1998). 

A. Wallraff et al, Nature 431, 162 (2004). 

L. Tian et al, Phys. Rev. Lett. 92, 247902 (2004); A. S. 

S0rensen et al, Phys. Rev. Lett. 92, 063601 (2004). 

P. Rabl et al, Phys. Rev. Lett. 97, 033003 (2006). 

K. Tordrup et al, Phys. Rev. A, 77, 020301 (2008); A. 

Imamoglu, Phys. Rev. Lett. 102, 083602 (2009); J. Verdii 

et al, Phys. Rev. Lett. 103, 043603 (2009); J. H. Wesen- 

berg et al, Phys. Rev. Lett. 103, 070502 (2009). 

D. I. Schuster et al, Phys. Rev. Lett. 105, 140501 (2010); 

Y. Kubo et al, Phys. Rev. Lett. 105, 140502 (2010). 

F. Jelezko and J. Wrachtrup, Phys. Stat. Sol. (a) 203, 
3207 (2006). 

J. M. Taylor et al. Nature Phys. 4, 810 (2008). 

J. E. Mooij et al. Science 285, 1036 (1999); C. H. van 

der Wal et al. Science 290, 773 (2000). 

See accompanying EPAPS Document. 

P. Bertet et al, Phys. Rev. Lett. 95, 257002 (2005); F. 

Yoshihara et al Phys. Rev. Lett. 97, 167001 (2006); K. 

Kakuyanagi, et al, Phys. Rev. Lett. 98, 047004 (2007); 

K. Semba, private communication (2009); M. Steffen et 

al, Phys. Rev. Lett. 105, 100502 (2010). 

H. K. Cummins et al, Phys. Rev. A 67, 042308 (2003). 

J. R. Maze et al, Phys. Rev. B 78, 094303 (2008). 

G. Balasubramanian et al. Nature Mater. 8, 383 (2009). 
M. D. Lukin, Rev. Mod. Phys. 75, 457 (2003); K. Ham- 
merer et al, Rev. Mod. Phys. 82, 1041 (2010). 

C. E. Lopez et al, Phys. Rev. A 75, 033818 (2007). 

H. Christ et al, Phys. Rev. B 75, 155324 (2007). 

J. M. Taylor et al, Phys. Rev. Lett. 90, 206803 (2003). 
F. C. Waldermann et al. Diamond and Related Materials 



5 



16, 1887 (2007). [23] J. Twamley and S. D. Barrett, Phys. Rev. B 81, 
[21] L. Childress et al, Science 314, 281 (2006); M. V. Gu- 241202(R) (2010). 

rudev Dutt et al, Science 316, 1312 (2007). 
[22] B. Smeltzer et al, Phys. Rev. A 80, 050302(R) (2009). 



Coupling Nitrogen Vacancy Centers in Diamond to Superconducting Flux Qubits: 

Supplementary Information 



D. Marcos,! M. Wubs,^'^ J. M. Taylor,'' R. Aguado,i M. D. Lukin,^ and A. S. S0rensen2 

' Departamento de Teoria de la Materia Condensada, 
Instituto de Ciencia de Materiales de Madrid, CSIC, Cantoblanco 28049, Madrid, Spain 
^ QUANTOP, The Niels Bohr Institute, University of Copenhagen, 2100 Copenhagen 0, Denmark 
'^DTU Fotonik, Danish Technical University, 2800 Kgs. Lyngby, Denmark 
Joint Quantum Institute and the National Institute of Standards and Technology, College Park, MD 20742 
''Department of Physics, Harvard University, Cambridge MA 02138 

, In the main text we have omitted several details from the estimates of decoherence and the 

performance of the procedures. In the following we give more details on how these were obtained. 
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I. DECOHERENCE OF THE NV CENTERS 



A. VSWAP operation between individual NV centers 

13 . 

• In the main text, we describe a method to perform an entangling operation between two NV centers near the same 
c/j ' FQ using this as a quantum bus. The interaction between different NVs is mediated off-resonantly by the FQ and 
D , can be described starting from the Hamiltonian HimA in the main text, extended to the case of 2 NVs. Taking the 
mixing angle 9 — > 7r/2, such that the FQ is optimally biased with respect to low-frequency flux noise, and a rotating 
frame such that the FQ is detuned —5, rather than the NVs being detuned 5, this Hamiltonian reads 



H = -'-if, + 1) + 'f. ( <?^W + <?a^)) + H.C. (1) 



^ ' In this picture, it becomes clear that there arc two zero-energy cigenstatcs: _ff|00)|0)i?Q = H\D)\Q) pQ ~ 0, with 

Q ; the dark state \D) = {w^^\lQ) - W^^|01)^ and the collective coupling constant for the two NVs G2 = 

: ,J{W[^)^ + (VFf '^)2/V2 = y^g'i + gl There is a corresponding bright state \B) = (w[^\lO) + Wf f |01)), 

' and the state |i?)|0)_FQ is coupled to |00)|1)_fq with coupling strength G2. Similarly the state |11)|0)fq is coupled to 

QQ \B)\1)fq with coupling strength G2. 

\ Therefore, in second order perturbation theory, projecting onto the lower energy state of the flux qubit, we have 
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H,ff = ^m{B\ + \n){n\). (2) 

' Here we will assume W^f « ^2 ~ V^g. The probability of finding the system in the excited state 

is then zero for the two dark states and 2g^/5^ for the two bright states of the system. Thus, FQ relaxation and 
Markovian dephasing, which occur directly to the FQ at a rate I/tJ'^, will lead to errors at a rate ^ 2g^/(5^Tj^ for 
, y ■ the two NV spins. Low- frequency noise enters in small variations of 5, which also contribute at order 2g^/S^. Since 
, such low-frequency noise is non-Markovian it will only give a quadratic contribution ~ [g'^t/ (S'^T^'^Yf' , which will be 
' less severe than the Markovian dephasing. We therefore take a worst-case scenario and assume the decoherence to be 
produced by Markovian dephasing. 

To optimize the gate operation, we note that the overall time when the FQ is coupled is tx = f 2p'' There is an 

additional wait time for the NV centers of 2tw ~ tx{^ — 40/7r), where <j) = 2 sin^^(l/-\/8) as in the main text. For 
NV spins in the states |00), \ D), the error is entirely from the total time tx + 2tvi/, as the FQ is not included, while 
for NV spins in the states \B), both FQ and NV dephasing errors enter. We take the latter case to overestimate 
the error induced in the operation, and find an effective dephasing of the NV center expf— (^^S^w^)'^], as expected 

from a dipolc-dipolc bath dccorrelation, and an additional dephasing due to the admixture of the FQ as cxp[— ^^^.fq ]■ 
The only free parameter we can optimize over is the detuning 5. Rewriting the total fidelity as 
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with a = ^(\^_20/7r) aiid P = the optimum occurs for 5^ = {a^P/2,f/^, giving i^, = cxp[-(/?/a)3/4(3i/4 + 3-3/4)j^ 

Substituting a and /3 in 5^ and the numbers given in the main text we find = cxp[— 2. 
and 5, = 0.96[/(T2NV)3/j.FQji/4 _ x 3.6 MHz. 

An additional error arises when Wl^^ ^ W['^, since then the unitary operation is not performing exactly the desired 

evolution; this leads to a reduced fidelity of the entangling operation by a factor W^fW^'^/g^. It is important to 
note, however, that this imperfection is not a decoherence effect and it may be possible to exploit the resulting unitary 
evolution even despite this imperfection, i.e., the resulting operation still resembles a vSWAP. 



B. Decoherence of the ensemble of NV centers 



In the main article we described the transfer of an excitation from the FQ to the collective state of an ensemble 
of NV centers. In this section we estimate the dephasing of this collective state of the ensemble due to magnetic 
dipole-dipole interactions with other spins in the diamond crystal. The paramagnetic impurities in the diamond 
crystal consist of both other NVs centers with either the same or different orientation as well as unpaired electronic 
spins on substitutional Nitrogens. In most experiments the ratio of Nitrogen to NV centers is typically quite low, e.g. 
on the order of 1% to 10% Below we refer to this ratio as the conversion efficiency. 

The dipole-dipole interaction is described by the Hamiltonian 0] 
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^dip 
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where rjk is the distance between the dipoles j and k, ejk is a unit vector between them, and m^^^ and m^'^^ are their 
magnetic moments. These can be expressed in terms of the spin operators of the impurities, but the quantization 
axis for the different spins will be different: NV centers will be aligned along the symmetry axis of the center whereas 
electron spins on Nitrogens will be quantized along the axis of the applied static field. To describe this situation we 



introduce spin operators si-'\ s^' and s);^' , defined relative to the quantization axis of each particular spin, and take 
the different orientations into account in the coupling constant describing the interaction between different spins. 

For the NV centers we furthermore ignore the ms = — 1 state which is assumed to be shifted out of resonance by a 
magnetic field. All the spins are therefore two-level systems which can be described by Pauli matrices and 
(7^'' . The spin operators can then be expressed as 



s(i) 
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(5) 



where the quantity Ij, which is unity for NV centers and vanishes for Nitrogen spin, accounts for the fact that the 
Nitrogen has spin 1/2 whereas the two-level NV system is made from the ms = and ms = 1 states of a spin-1 
particle. The magnetic dipole Hamilonian is then given by 



dip 
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(6) 



where we have used the rotating- wave approximation to ignore terms which do not conserve energy. By this approx- 
imation, the first term vanishes between spins with different resonance frequencies, i.e., with different orientations. 
Notice, that if the quantization axes for spin j and k are different, then the angle between these axes can be contained 
in the coupling constants ajk and bjk- This is the case for the interaction between the NV and Nitrogen spins, whose 
quantization axes, determined by the crystal axis and the external magnetic field respectively, form an angle /3 between 
them. For the sake of simplicity, we here only present the case f3 = 0. This is a worst-case scenario for the Tj* derived 
below, but the general situation /? ^ only leads to minor modifications. With parallel quantization axes, we have 
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where Qjk is the angle between the vector Cjk and the crystal axis, e^'j^ := (gJ'^ ±ief )/2, and 
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with the Bohr magneton and ge"'^ the electron g-factor of the jth spin. 

To evaluate the dephasing of the spin wave we will assume that we start out in a state |00...0)nv where all the NVs 
are initially prepared in their ground states, and that a superposition state co|0)fq + ci|1)fq is transferred into the 



collective spin wave, resulting in a state (co + ci J+)|00...0)nv, with J+ = (l/G) J2j 
evolution of the coherence 



(i) 



We then evaluate the time 



(J_(0> = Tr[e*^*J_e-*^* 



(co + ci J+)|00...0)nv(0...00|(cS + J_)pe 



(9) 



where pb is the initial density matrix of the bath. Since the Hamiltonian ^ conserves the number of excitations in 
the NVs, this expression can be simplified to 



(J_(t)) = cSciTr[e*^*J_e-*^*J+|00...0)Nv(0...00|pB]. 



(10) 



Writing out the above equation in terms of single-spin operators we see that the dephasing is determined by the 

time evolution of the two-point correlation function (o'i^'' (t)(T^'' (i = 0)) in the state |00...0)nv(0..-00|/9b- Since the 
density of Nitrogen is much higher than the density of NVs, the dephasing of the spin wave will predominantly be due 
to the interaction with the electronic spins of Nitrogen atoms, and we shall therefore ignore the interaction among 
the NVs. In this approximation there is no longer a mechanism in the Hamiltonian which can transfer the excitation 
from one NV to another and the correlation function (ct^-' {t)a'^'^ (t = 0)) vanishes exactly for j ^ k. The dephasing 
of the spin waves thus reduces to the calculation of the single-spin dephasing averaged over the spin wave: 
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where the sum runs over all NVs in the spin wave. Assuming all NVs to be equivalent this can be simplified to the 
calculation of the dephasing of a single-spin in a bath 



(J_(t)) = (J-(t = 0))Tr 
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When considering a single NV, the action of the spin bath described by ([6|) essentially corresponds to a random 
magnetic field generated from the spin ensemble along the crystal axis. This effective field will be fluctuating in time 
because the dipole-dipole interaction among the Nitrogen spins introduces flip-flop processes. When the Nitrogen 
concentration is much higher than that of NVs, the distance between Nitrogens is comparable to the distance between 
NVs and the nearest Nitrogen. These flip-flop processes will therefore take place on a time scale which is comparable 
or only slightly faster than the NV dephasing time. For simplicity we here ignore the flip-flop processes and consider 
a static environment. This represents a worst-case scenario, since this approximation removes the time averaging of 
the fleld from these processes. The actual dephasing time Tj* will therefore be slightly larger than what we predict 
here. Within this approximation the expression in Eq. (|12[) can be reduced considerably: 



(J_(0) = (J-(t = 0))Tr 
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(13) 



where the product is over the Nitrogen spins. 

With expression (jl2p we can evaluate the dephasing for a given bath. For a given density n, the typical strength of 
the spin-spin interaction is fiQ^'^g'^n/ATT, of the order of ^ 25 for high Nitrogen densities tt-n ~ 10^^ cm~'^. This 
is much smaller than the typical operating temperature of the FQ (tens of mK), and we can therefore neglect the 
interaction among the spins for determining the initial density matrix of the bath pB ■ Furthermore any no n- vanishing 
mean value (o'^'''') 7^ only results in a mean shift of the resonance frequency of the ensemble. The dephasing will 
thus be determined by spin fluctuations. A worst-case scenario can be obtained by assuming that the mean value 



vanishes (a 



0, in which case the variance is maximal. The evolution of the coherence is then given by 



(J_(i)) = (J_(i = 0))Y[cos{4b,kt). 
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To simplify this expression we expand it in time and find 
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where we have replaced the sum by an integral and introduced the Nitrogen density n-^. Using Eq. ([7]) this expression 
reduces to 



(J_(0) = (J-(< = 0)> 
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The angular integral gives 8/5 but the radial integral has a strong divergence at ?■ — > 0. Since the integral is an integral 
over the distance between the NV and the Nitrogen impurity, this divergency represents the very fast dephasing of NV 
centers which happen to have a nearby Nitrogen spin. Such NV centers with a nearby Nitrogen spin will, however, 
also be far out of resonance during the interaction with the Flux qubit. This interaction takes a time l/G, and 
therefore any NV with a dipole interaction stronger than G will effectively not participate in the spin wave. We 
exclude these NVs by truncating the integral at a distance rmin when the interaction strength reaches the value of 
the coupling constant G, i.e., 

^"84^^^3 ■ ^^^> 
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This truncation of the integral gives 

(J_(t)) = (J_(f = 0)) 
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where rj is the Nitrogen-to-A^V^ conversion efficiency. The result corresponds to the first term in an expansion of 
an exponential decay exp(— (i/Tj*)^) with coherence time 

T; = . (19) 

V 15^ "NV<^ 

Notice that in order to be applicable, our regularization procedure requires rmin to be much smaller than the typical 
distance between spins 1/ -^/nN, since we require that only a small fraction of the NVs are excluded. This condition 
is fulfilled in the interesting regime GTj* ^ 1. Taking a conversion efficiency 77 = 0.05, we find Tj* « 0.3 iJS for a 
FQ with L = 5 /xm and a density unv = 10^^ cm"^, corresponding to G « 27r x 15 MHz. For a full transfer of the 
state from the FQ to the spin wave and back we need a transfer time t = tt/G, corresponding to an infidelity of the 
order of 1 — ~ [1/2X2)'^ < 0.5% (the factor of 2 accounting for the fact that the excitation only spends half of the 
transfer time in the spin wave). For different densities the infidelity scales as 1 — cx ^/n^v, becoming smaller at 
lower densities due to the reduced dipole-dipole interaction, e.g., at riNv = 10^^ cm~^ the error is reduced by an order 
of magnitude but then we are approaching the limit where the FQ decoherence becomes important. The transfer of 
excitations from the FQ to the spin wave is thus feasible in the regime 10^^ cm~'^ < tt-nv ^ 10^* cm~'^. In order for 
the spin system to be useful as a long-term memory the coherence time should, however, also be sufficiently long to 
allow for the transfer to the nuclear spin for long-term storage. Since this can at best be achieved on a time scale set 
by the hyperfine interaction (^ 5 MHz), this would exclude working at the highest densities in this interval. Working 
at a density of unv ~ 10^^ cm~'^ leads to Tj* w 1.8/iS, which is sufficient to allow a transfer of the excitation from the 
electron spin to the nuclear spin with an infidelity 1 — at the percent level. 

The estimates above indicate that it is realistic to achieve a transfer of excitations from the FQ to the spin wave 
and back even without extending the coherence time by spin-echo techniques. These may, however, be desirable in 
order to achieve even longer coherence times. In particular T2 may be extended if the NVs can be flipped by an 
external AC field on a time scale faster than the Nitrogen flip-flop processes. One should, however, be careful about 
applying spin-echo to the collective NV spin since, e.g. errors in the pulse area may give rise to collective decoherence 
processes. These can be more detrimental to quantum states stored in collective degrees of freedom than in individual 
spins (for instance flipping the |0) — |1) transition would require control of the pulse area to an accuracy better than 
1/VN in order to preserve the collective state). Since the dephasing of the NV spin wave is dominated by Nitrogen 
spins a more desirable solution could be to apply the external driving field to the electron spins on the Nitrogen 
atoms. If these spins are flipped on a time scale much faster than the flip-flop processes it would lead to an increased 
coherence time of the NV spin wave. Since Nitrogen spins have a considerably different resonance frequency, this 
could be achieved with little influence on the NV spin wave where the quantum information is stored. 



II. DECOHERENCE OF THE FLUX QUBIT DUE THE SPIN BATH 



A different concern is the extent to which the diamond crystal with a high density of spins may cause dephasing of the 
flux qubit, of particular importance for the coupling with an ensemble of NVs. As discussed above, the paramagnetic 
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impurities in the diamond crystal consist primarily of NV centers and unpaired electron spins on Nitrogen atoms which 
were not converted into NVs centers during the annealing process. Due to the typically low conversion efficiency from 
Nitrogen to NV centers, we shall first consider the effect of Nitrogen impurity spins and then discuss the role of the 
NV centers. The coupling of the FQ to paramagnetic impurities can be described by a Hamiltonian similar to Eq. 
(1) in the main text: 



H = iJspin + eh/2 + Afi +hY. • "^^'^ 



(20) 



where the spin Hamiltonian Hapin = J2j ^jS^^ /2 

the interaction between them, encapsulated in iJjnt 
Hamiltonian is transformed into 



Hint describes an energy splitting Aj of the individual spins and 
Changing to the dressed-state picture of the flux qubit, this 



H = H. 
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In the rotating frame with respect to uf'^/2, the operator 7^3 remains stationary whereas fi oscillates at a frequency uj 
- much higher than any time scale of the bath, provided that there are no near resonance impurities in the diamond 
sample. This is the case for the Nitrogen spins, whose splitting is determined by the applied field. The NVs are near 
resonance and will be dealt with below. The slowly varying and rapidly oscillating terms have different qualitative 
behavior so wc will consider them separately. 

Let us first consider the slowly-varying contribution. As we will show now this contribution vanishes if we work 
close to the degeneracy point of the FQ cos k, Q where the left- and right-circulating current states are degenerate. 
The slowly- varying contribution is described by 



low 



(22) 



with 
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and I3j being the angle between the quantization axis the jth spin and the crystal axis of the NV (recall that the 
z-axis is defined by the crystal axis). For simplicity, we define 



W-f «(r,) cos(/3,) - VKf «(r,) sin(/3,). 



(24) 



Since the interaction of the flux qubit with an individual spin is on the order of 5, the former has little influence on 
the state of an individual spin for the duration of the interaction ^ 1/G ^ 1 / g^/N . Furthermore the coupling of the 
FQ to collective degrees of freedom in the spin bath, will only have a limited influence on the state of an unpolarized 
bath, since the FQ can fiip at most a single-spin. We can therefore ignore the influence of the FQ on the spin bath 
and consider W^^^ as a fluctuating external field. Any mean value of this field will merely give rise to a shift of the 
resonance frequency which can be compensated and we thus need to consider the fluctuations. The root mean square 
of the fluctuating field is given by 



5W^^ = cos 



1/2 



(25) 



where we have used the cumulant notation {{ab)) :~ {{a — {a)){b — (6))). As discussed in the previous section, the 
temperature is typically high compared to the dipole-dipole interaction energy. We can therefore ignore correlations 
of different spins and consider only the j = k contribution, giving 



J^'S ^ cos(0)[5:.|(((5l^') 



1/2 



(26) 



Considering a Nitrogen-to-NV conversion oirj ^ 5%, the sum over the impurities in the ensemble will be dominated by 
the Nitrogen impurities. An estimate of this expression can be obtained by noting that the coupling to the Nitrogen 
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spin K is comparable to the coupling g to the NVs. The difference between this sum and the sum leading to the 
collective coupling constant G only differ by the number of terms in the sum and we obtain 

^ ^-^G. (27) 

A more accurate treatment taking into account the full spatial distribution of the field over a sample of the dimensions 
considered in the text only changes the estimate by a factor of less than 2. Eq. ([27|) quantifies how close to the 
degeneracy point we need to be so that the dephasing of the FQ induced by the bath of spins is negligible. In 
particular, from pT]) we derive that we can safely neglect this dephasing on a time scale 1/G for a conversion 
efficiency tj = 0.05 if cos(6l) < 0.01, for which SW^^^/G < 0.05. 

Next we turn to the rapidly oscillating part of the coupling, described by the term containing fi, and assume 
that we are near the degeneracy point, that is sin0 « 1. We will now argue that since fi is oscillating rapidly, the 
slowly- varying dynamics of the spin bath will only have a very weak influence on the FQ. Specifically the spin bath can 
influence the FQ cither through direct transitions to the spin bath or through a slow dephasing. The direct coupling 
can be excluded by noting that the coupling constant to collective excitations of the bath is limited by G/y/fj, but for 
reasonably low applied magnetic fields, i?'^''' < 10 mT, the detuning oj — Aj will be of the order of the NV zero-field 
splitting (a few GHz). The probability to transfer the excitation is ^ G'^ /{lu — Aj)^?/, and can thus safely be neglected 
since G ~ MHz. The dephasing of the FQ caused by the bath can be estimated by calculating the energy shift of the 
FQ in second order perturbation theory. Assuming that apart from the free precession, the spins change slowly on a 
time scale set by w — Aj , the effective interaction is given by 

^eff = rsX^ 4^ + 4(^2 _A2) - 2(c.2-A2) ^^^^ ^^^^ 

Here the first two terms are independent of the state of the impurities, and therefore merely lead to a mean shift of 
the energy, which can be compensated by a magnetic field. It is thus only the last term which leads to dephasing. 
Again, we can estimate the root mean square value of this term by neglecting the correlations among the impurities 
and we find a typical energy shift 

which we can safely neglect since the single-spin coupling g ^ kHz is much smaller than any other quantities in the 
system. 

A remaining problem is the influence of the NVs with different orientations. The interaction with these centers 
is similar to that with the resonant subensemble and of order G. We assume that we select a single orientation 
by applying a magnetic field with a component along the axis of the NV center. This will shift all other centers 
out of resonance by an amount Su ~ gefisB^^^, such that the resulting error can be estimated to be G^/Suj'^. For a 
reasonable applied field, e.g., B°^*^ ^ 10 mT, this error is negligible. 
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